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Abstract: The whole mitochondrial genome sequence of red fox (Vuples vuples) was determined. It had a total 


length of 16 723 bp. As in most mammal mitochondrial genome, it contained 13 protein coding genes, two ribosome RNA 
genes, 22 transfer RNA genes and one control region. The base composition was 31.3% A, 26.1% C, 14.8% G and 27.8% 
T, respectively. The codon usage of red fox, arctic fox, gray wolf, domestic dog and coyote followed the same pattern 
except for an unusual ATT start codon, which initiates the NADH dehydrogenase subunit 3 gene in the red fox. A long 
tandem repeat rich in AC was found between conserved sequence block 1 and 2 in the control region. In order to confirm 


the phylogenetic relationships of red fox to other canids, phylogenetic trees were reconstructed by neighbor-joining and 


maximum parsimony methods using 12 concatenated heavy-strand protein-coding genes. The result indicated that arctic 
fox was the sister group of red fox and they both belong to the red fox-like clade in family Canidae, while gray wolf, 
domestic dog and coyote belong to wolf-like clade. The result was in accordance with existing phylogenetic results. 
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Mammalian mitochondrial genomes typically have 
a set of 13 protein-coding genes, two ribosome RNA 
genes (12S RNA and 16S RNA) and 22 transfer RNA 
genes. The gene order is highly conserved among most 
vertebrates (Boore, 1999). The two strands that make up 
the genome are commonly known as the heavy strand 
(H-strand) and the light strand (L-strand). Twelve of the 


Received date: 2009-09-11; Accepted date: 2010-01-20 





XERRA: A 


MBWS: 0254-5853-(2010)02-0122-09 


13 protein-coding genes locate on the H-strand and only 
ND6 gene are located on the L-strand. Vertebrates 
usually contain two non-coding regions, the major of 
which 
transcription promoters. Consequently, it is known as the 
control region (CR) (Wolstenholme, 1992; Boore, 1999). 
The small non-coding region is the origin of replication 


contains mitochondrial replication and 
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of the L-strand (OL), thought to have a functional role in 
replication (Shadel & Clayton, 1997). 

Belonging to Canivora, Canidae, Vuples, red fox 
(Vuples vuples) is one of the worldwide distributed 
mammals. Comparing complete animal mitochondrial 
genome sequences is now common for phylogenetic 
reconstruction and as a model for genome evolution 
(Wei et al, 2008). Mitochondrial DNA sequences were 
also used to clarify phylogenetic relationships within 
Canidae (Geffen et al, 1992; Wayne et al, 1997). 
Sequence data from three mitochondrial genes (Wayne et 
al, 1997) suggest that Canidae falls to four monophyletic 
groups: (1) the wolf- and jackal-like canids; (2) the red 
fox-like canids; (3) the South American foxes; and (4) 
the maned wolf (Chysocyon brachyurus) and bush dog 
(Speothos venaticus). Analysis of various morphological 
and mitochondrial DNA data (Zrzavy & Ricankova, 
2004) as well as phylogenetic analysis using six nuclear 
loci combined with mitochondrial data (Bardeleben, 
2005) both agree with the three clades, including the red 
fox-like canids, the South American foxes, and the 


wolf-like canids. 

This paper reports a complete mitochondrial 
genome of red fox, which was also compared with that of 
other canids to discuss the red fox mitochondrial genome 
structure and evolution. A phylogenetic tree was 
reconstructed to confirm the evolution position of red fox 
in Canidae. 


1 Materials and Methods 


1.1 Samples and DNA extracting 

A blood sample of 200 uL from a female red fox 
was obtained from Beijing Zoo. Total genomic DNA was 
extracted following the method of Sambrook & Russell 
(2001) and the obtained DNA solution was preserved 
under —20 °C for next use. 
1.2 PCR and sequencing 

Five pairs of primers were designed based on the 
sequences of wolf and coyote, which were available 
online with the accession number of NC008092 and 
NC008093, respectively. Primer sequences are shown in 
Tab. 1. 


Tab. 1 Five primer pairs for the first PCRs 


Primer pair Primer Sequence(5'—3') 
1 1-F TCCCTCTAGAGGAGCCTGTTC 
1-R GGGTATGGGCCCGATAGCTT 
2 2-F GGCGGATAAAAGAGTTACTTTGATAGAG 
2-R GCGAATTTAACTTTGACAAAGTCATGT 
3 3-F GAAGAAAGGAAGGAATCGAACC 
3-R GCGTAGGGATGATAATTTTTAGCATT 
4 4-F GTATTTGCTGCCTGCGAAGC 
4-R TAGTGGTGGGATTGGTTGTGC 
5 5-F GGGTATTGCTCAGTAGCCATAGC 
5-R GGTTTGCTGAAGATGGCGGTATAT 


PCR was implemented in a 50 uL system including 
1 uL DNA, 5 uL 10 x Long-PCR Buffer, 4 uL dNTP 
mixtures (2.5 mmol/L each NTP), 15 mmol MgCl; pfu 
polymerase 1 uL (5 U/L); 2 uL primers (10 pmol/L). 
After degeneration of 95°C for 5 min, 35 cycles of 95°C 
degeneration for 30 s, 57°C annealing for 30 s, and 72°C 
extension for 4.5 min were ran. Final extension was ran 
at 72°C for 10 min. PCR products were purified from a 
1% agarose gel and were then used as DNA templates in 
secondary PCR. Secondary PCR was conducted with 10 
pairs of primers (Tab. 2). PCR conditions were: initial 
degeneration of 95°C for 5 min, then 35 cycles of 95°C 
degeneration for 30 s, 57°C annealing for 30 s, 72°C 
extension for 2.5 min and a final extension step of 72°C 


for 4.5 min. Reacting system and purifying methods are 
the same as the first PCR. The purified PCR products 
were then sent out for sequencing. 
1.3 Sequence assembly, annotation and analysis 
Three times of sequencing were performed to make 
sure correct target sequences were amplified. The 
DNAMAN 6.0 program (Lynnon Biosoft, Quebec, 
Canada) was used for editing and assembling the 
sequence. The positions of protein coding genes, rRNA 
genes and non-coding regions were located through 
alignment with that of wolf and coyote using BLAST 
(http://blast.ncbi.nlm.nih.gov/Blast.cgi). Start and stop 
codons were identified using the vertebrate mitochon- 
drial code except some potential incomplete stop codons. 
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Tab. 2 


Primer pair Primer 


Ten primer pairs used for secondary PCRs 


Sequence(5'—3') 





1 1-F 


TCCCTCTAGAGGAGCCTGTTC 


1-R TCCGAGGTCACCCCAACC 


2 2-F 


GACGAGAAGACCCTATGGAGC 


2-R GGGTATGGGCCCGATAGCTT 


3 3-F 


GGCGGATAAAAGAGTTACTTTGATAGAG 


3-R GCCTACTATACCGGCTCATGC 


GCTCAGCCATTTTACCTATGTTC 


4-R GCGAATTTAACTTTGACAAAGTCATGT 


5 5-F 


GAAGAAAGGAAGGAATCGAACC 


5-R GCGAAGAGTTGTAGTGAAATCATAT 


GCTACCTAATGACCCACCAAAC 


6-R GCGTAGGGATGATAATTTTTAGCATT 


7 7-F 


GTATTTGCTGCCTGCGAAGC 


7-R CGCTTATCTGGAGTTGCACC 


8 8-F 


CCGCAAGAACTGCTAATTCATG 


8-R TAGTGGTGGGATTGGTTGTGC 


9 9-F 


GGGTATTGCTCAGTAGCCATAGC 


9-R GCAGAATTTCAGCTTTGGGTG 


10 10-F 


CGCGATGAAGAGTCTTTGTAGTAT 


10-R GGTTTGCTGAAGATGGCGGTATAT 


21 tRNAs genes were identified by tRNAscan-SE Search 
Server v.1.21 (http://lowelab.ucsc.edu/tRNAscan-SE/). 
The tRNASer (AGC) gene not found by tRNAscan-SE 
was determined by comparison with the two canids 
mentioned above. RNA structure 4.6 (Lowe & Eddy, 
1997) was used to speculate the second structure of the 
22 tRNAs. MEGA v.4.0 (Tamura et al, 2007) was used to 
compute base composition, substitution and codon usage. 

The mitochondrial genome sequence of the red fox 
was submitted into GenBank with the accession number 
of GQ374180. 
1.4 Phylogenetic analysis 

Four complete protein sequences available for 
Canidae were downloaded from GenBank, and tiger 


(Panthera tigris) was used as an outgroup for 
phylogenetic analyses (Tab. 3). All protein coding genes 
except the ND6 gene were chosen to reconstruct 
phylogenetic tree. Multiple alignments were performed 
by ClustalX 1.83 (Thompson et al, 1997) with the default 
parameters. The results of alignments were also carefully 
checked and edited by eye. Phylogenetic analyses were 
conducted using neighbor-joining (NJ) and maximum 
parsimony (MP). The NJ analysis was performed in 
MEGA 4.0 using Kimura 2-parameter as the nucleotide 
substitution model. Statistical confidence was assessed 
by bootstrap analysis with 1000 replications. The MP 
analysis was performed in PAUP*4.0b10 (Swofford, 
2003) using a heuristic search. The robustness of the 


Tab.3 The taxa and sequence’s accession number of the six species used for phylogenetic analyses 





Species Common name Accession No. and references 

Vuples vuples Red fox GQ374180,This paper 

Alopex lagopus Arctic fox AHO014073, Delisle & Strobeck, 2005 
Canis lupus Gray wolf NC_008092, Bjornerfeldt et al, 2006 
Canis lupus familiaris Dog NC_ 002008, Kim et al, 1998 

Canis latrans Coyote NC_008093, Bjornerfeldt et al, 2006 
Panthera tigris Tiger EF551003, unpublished 
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tree was tested with 100 bootstrap replications. 


2 Results 


2.1 Mitochondrial genome structure 

The total length of red fox mitochondrial genome is 
16 723 bp, a little shorter than that of gray wolf(1 6729 
bp) and coyote (16 724 bp). The red fox mitochondrial 
genome shares high similarity with those of three other 
canids: gray wolf (84.93%), domestic dog (84.86%) and 
coyote (84.87%). Similar to other animals (Boore, 1999), 
the red fox mitochondrial genome has 13 protein-coding 
genes, 22 tRNA genes, two rRNA genes (12S rRNA, 16S 
rRNA) and one control region (CR). Except for the 
protein coding gene ND6 and 8 tRNA genes [tRNAGIn, 
tRNAAla, tRNAAsn, tRNACys, tRNATyr, tRNASer 
(UCN), tRNAGlu, tRNAPro] encoded at the L-strand, 
the rest genes are encoded on the H-strand. Overall 
organization of the genome is comparable to the situation 
of other vertebrates(Tab. 4). As shown in Fig. 1, the 
genes are so compactly arranged that few gaps were 
found between them and some genes even overlap each 
other. The longest gap was found between COII and 
tRNALys by 16 nucleotides, while ATPase8 and 
ATPase6 overlap each other by 33 nucleotides, which is 
the longest one. 

The overall base composition of the red fox 
mitochondrial genome is 31.3% A, 26.1% C, 14.8% G 
and 27.8% T, respectively (Tab. 5).The consequent order 
would be A>T>C>G, which was in coincidence with the 
order of other mammals. The A+T content is 59.1%, 
which indicated an A-T rich in the red fox mitochondrial 
genome. 

2.2 Protein coding genes 

Except ATA for ND2, NDS5 and ATT for ND3, other 
protein coding genes use the start codon ATG. When it 
comes to stop codons, NDI, COI, COIL, ATPase8, 
ATPase6, ND4L, ND5 and ND6 terminate with TAA, 
while ND2 and Cytb terminated with TAG and AGA 
respectively. Besides, three incomplete stop codons were 
detected: T (COM), TA (ND3) and T (ND4). Codon 
usage analysis (Tab. 6) indicates obvious codon bias and 
CUA (Leu), AUA (Met), AUU (lle) are the most used 
three codons. Codon degenerate was also found by codon 
usage analysis. Both four-fold degenerate codons (e.g. 
codons specify Ala and Pro) and two-fold degenerate 
codons (e.g. codons specify Asn and Tyr) were found. 
Leu and Ser are found to be specified by both four-fold 
degenerate codons and two-fold degenerate codons. 

The total length of 13 protein coding genes is 


11 411 bp, accounting for 68.24% of the complete 
sequence of the red fox mitochondrial genome. As shown 
in Fig.1, three pairs of protein coding genes, ATP8/ATP6, 
ND4L/ND4 and NDS5/ND6 overlap each other by 33 bp, 
7 bp and 17 bp respectively. For the 13 protein-coding 
genes of the red fox mitochondrial DNA, the nucleotide 
distribution at three codon positions differ significantly 
(P<0.01).A strong bias against G (only 9.2%) at the third 
codon position was observed. The frequency of A+T in 
all protein-coding genes ranges from 53.6% to 61.6%. 
2.3 Ribosomal RNA and transfer RNA genes 

The 12S rRNA gene of the red fox mitochondrial 
genome was located between tRNA™ gene and tRNA“ 
gene with a length of 957 bp, longer than that of gray 
wolf (954 bp) and coyote (955 bp) due to different 
numbers of insertions. The base composition is 35.6% A, 
22.2% C, 18.6% G and 23.6% T respectively. The 16S 
rRNA gene was located between tRNA” gene and 
tRNA‘) gene. The length is 1579 bp, which was 
identical to coyote but shorter than gray wolf (1580 bp). 
The base composition is 36.8% A, 20.6% C, 17.6% G 
and 25.0% T respectively. The computed A+T content of 
the 12S rRNA and 16S rRNA gene are 59.2% and 61.8%, 
which also indicates an A-T rich. 


Vuples vuples 


16,723 bp 





Fig. 1 Organizations of red fox mitochondrial 
genome. Transfer RNA genes are depicted 
by their one-letter amino acid codes 
Numbers indicate non-coding nucleotides between genes or gene 
overlap (negative values). Arrows indicate orientation on (+) strand 


(clockwise) or (—) strand (counterclockwise). 
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Tab.4 Organizations and characteristics of red fox mitochondrial genome 











Position number 5 Codon ; 
Name of gene ———————_____ Size(bp) _ Intervals Anti-codon Strand 
Start Stop Start Stop 

tRNA 1 69 69 0 GAA H 
12SrRNA 70 1026 957 0 H 
tRNA 1027 1093 67 0 TAC H 
16S rRNA 1094 2672 1579 0 H 
tRNA" UUR 2673 2747 75 2 TAA H 
ND1 2750 3706 957 0 ATG TAA H 
tRNA"? 3706 3774 69 -3 GAT H 
tRNA®™ 3772 3845 74 1 TTG L 
tRNA™* 3847 3916 70 0 CAT H 
ND2 3917 4960 1044 -2 ATA TAG H 
tRNA™ 4959 5026 68 12 TCA H 
tRNAA® 5039 5107 69 1 TGC L 
tRNAA" 5109 5179 71 0 GTT L 
OL 5180 5217 38 -4 L 
tRNA“S 5214 5281 68 0 GCA L 
tRNA™ 5282 5349 68 1 GTA L 
COI 5351 6895 1545 -3 ATG TAA H 
tRNAS*UN) 6893 6961 69 6 TGA L 
tRNA“? 6968 7035 68 0 GTC H 
COII 7036 7719 684 16 ATG TAA H 
tRNAs 7737 7803 67 1 TTT H 
ATPase8 7805 8008 204 -33 ATG TAA H 
ATPase6 7966 8646 681 -1 ATG TAA H 
COII 8646 9429 0 ATG T H 
tRNASY 9430 9497 68 0 TCC H 
ND3 9498 9844 347 0 ATT TA H 
tRNA“® 9845 9913 69 0 TCG H 
ND4L 9914 10210 297 -7 ATG TAA H 
ND4 10204 11582 1379 -1 ATG T H 
tRNA 11582 11650 69 1 GTG H 
tRNA SAY) 11652 11710 59 0 H 
tRNAEUCON) 1711 11780 70 0 TAG H 
NDS 11781 13601 1821 -17 ATA TAA H 
ND6 3585 14112 528 0 ATG TAA L 
tRNAS" 14113 14181 69 4 TTC L 
Cytb 4186 15325 1140 0 ATG AGA H 
tRNA™ 15326 15395 70 0 TGT H 
tRNA?” 5395 15460 66 0 TGG L 
D-loop 5461 16723 1263 0 


The length of the red fox mitochondrial tRNA genes stem and loop (DHU Stem and loop), other tRNAs can 
range from 59 bp to 75 bp. Among the adjacent tRNA fold into typical cloverleaf secondary structure (Fig.2). 
genes, only tRNAIle and tRNAGIn overlap each other by The average A+T content in all tRNA genes is 63.6%, 
3 bp. Except for tRNASer(AGY) lack the dihydrouridine higher than that of protein-coding and rRNA genes (Tab. 
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Fig.2 Typical cloverleaf secondary structure (Left, e.g. tRNA’) and the secondary structure of RNAS 4°”) (Right) 
Six major elements are included in the typical structure of tRNA: Anticodon stem and loop, T¥C Stem and loop, Amino Acid Acceptor 


Arm, Extra Loop, Discriminator Nucleotide and finally a DHU Stem and loop. Whereas the secondary structure of tRNA‘ “® lack of the 


DHU Stem and loop. 


Tab. 5 Comparison of base compositions in different genes (or regions) of the red fox mitochondrial genome 


Gene/region 

A C 
Protein coding 
Ist 31.1 24.7 
2nd 19.4 26.2 
3rd 36.6 30.1 
Total 29 24 
tRNA 34.7 20.6 
12S rRNA 34.7 22.2 
16S rRNA 34.7 20.6 
D-loop 34.7 29.2 


Base composition (%) 


G T A+T 
21.8 22.5 53.6 
12.2 42.2 61.6 
9.2 24.1 60.7 
14.4 29.6 58.6 
15.9 28.9 63.6 
18.6 23.6 59.2 
17.6 25 61.8 
15.8 25.1 55.5 
14.7 27.7 59.1 


Overall 34.7 26.1 


5). 
2.4 Non-coding regions 

Locating between tRNA”? and tRNA’ genes, the 
CR of the red fox mitochondrial genome has a length of 
1263 bp, shorter than that of gray wolf (1269 bp) and 
coyote (1266 bp). The base composition is 29.9% A, 
29.2% C, 15.8% G and 25.1% T, respectively. The 
sequence (ETAS) 
located at the 5' end of the CR and contain a termination 


extended termination associated 


associated sequence (TAS). This domain contains a core 
motif of and its complement and fold into a ‘hairpin’ 
structure; the central conserved domain (CD) domain 
includes five conserved blocks (CSB-F, E, D, C, B); the 
conserved sequence block (CSB) domain contains three 
conserved blocks (CSB-1, 2, 3). A region of tandem 
repeat region called RS3 (Douzery & Randi, 1997) was 
determined by Tandem Repeats Finder program (Benson, 
1999) at region 645—981bp. Similar to some mammals 
investigated so far (Fumagalli et al, 1996; Hoelzel et al, 


1994; Ketmaier & Bernardini, 2005), it was located 
between CSB-1 and CSB-2. Particularly, the base 
composition of A+C is 70.0%, indicating an A-C rich in 
this region. The single repeat is a motif of TACACACG 
(8 bp). There are totally 37 copies, 32 of which are 
complete repeats. Besides, two repeats show the same 
transition of AG, CT and G—A at position 4, 5 and 
8 respectively and the sequence is TACGTACA. Two 
repeats show the same transition of T—C, A—G and 
C—>T at position 1, 4, 5 respectively and the sequence is 
CACGTACG. One repeat shows a transition of A—>G 
and CT at position 4 and 5. Among this tandem repeat 
region, we found an alternative array of repeats. In this 
case there were 14 repeats with the motif of 
ACGTACACACGTACACACGT (20 bp). Except for 
one repeat with a transition of G>A, T—C, A—>G and 
CT at position 3, 4, 7 and 8 respectively, the other 13 
repeats follow the complete repeat. 
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Tab. 6 Codon usage and codon bias (RSCU) in red fox mitochondrial genome 
AA Codon Count RSCU AA Codon Count RSCU AA Codon Count RSCU AA Codon Count RSCU 
Ala GCA 97 1.45 Glu GAA 64 136 Lys AAA 84 165 Ter AGA 1 0.4 
GCC 93 1.39 GAG 30 0.64 AAG 18 0.35 AGG 0 0 
GCG 15 0.22 Gly GGA 82 153 Met ATA 189 1.51 TAA 8 3.2 
GCT 63 0.94 GGC 50 0.93 ATG 61 0.49 TAG 1 0.4 
Arg CGA 40 3.64 GGG 41 0.76 Phe TTC 137 1.2 Thr ACA 108 1.44 
CGC 7 0.64 GGT 42 0.78 TTT 92 0.8 ACC 101 1.35 
CGG 8 0.73 His CAC 75 1.46 Pro CCA 68 1.4 ACG 11 0.15 
CGT 10 0.91 CAT 28 0.54 CCC 68 1.4 ACT 79 1.06 
Asn AAC 84 1.16 Ile ATC 154 0.91 CCG 7 0.14 Trp TGA 88 1.69 
AAT 61 0.84 ATT 183 1.09 CCT 51 1.05 TGG 16 0.31 
Asp GAC 48 1.28 Leu CUA 230 2.34 Ser AGC 33 0.69 Tyr TAC 73 1.08 
GAT 27 0.72 CUC 94 0.96 AGT 21 0.44 TAT 62 0.92 
Cys TGC 19 1.52 CUG 63 0.64 TCA 86 1.8 Val GTA 74 1.58 
TGT 6 0.48 CUU 81 0.83 TCC 67 1.4 GTC 48 1.03 
Gin CAA 73 1.64 TTA 94 0.96 TCG 14 0.29 GTG 25 0.53 
CAG 16 0.36 TTG 27 0.28 TCT 66 1.38 GTT 40 0.86 


The putative light strand replication (OL) was 
detected in the region WANCY, which contains five 
tRNA genes: tRNA”, tRNA“", tRNA“, tRNA“ and 
tRNA™. This region is 38 bp long and overlaps the 
tRNA“ gene by 4 bp. It can fold into a stem-loop 
secondary structure with a 12 bp stem and a loop of 14 
nucleotides. 

2.5 Phylogenetic analysis 


100 


m 
0.02 


The same topology was obtained from both NJ tree 
and MP tree (Fig. 3). The bootstrap values are very high 
and all up to 100. The red fox and arctic fox fall into one 
clade, while gray wolf, domestic dog and coyote fall into 
another one. The red fox is a sister group with arctic fox. 
The coyote diverged after the emergence of red fox and 
arctic fox and prior to gray wolf and dog. 


100 Gray wolf 
100 Dog 
Coyote 
Red fox 
Arctic fox 


Tiger 


Fig.3 Phylogenetic relationships by neighbor-joining (NJ) and maximum parsimony (MP) methods inferred 


from concatenated 12 protein-coding genes 


Same topology was obtained and all bootstrap values are 100. 


3 Discussion 


3.1 Mitochondrial genome evolution 
The size and structure of mitochondrial genome of 
vertebrates are highly conserved. The length ranges from 


15.7 to 19.5 kb and the red fox, gray wolf and coyote all 
fall within this. There are only several nucleotides 
differences in length among these three canids. The 
arrangement of multiple genes in red fox is in line with 
most vertebrates and the overall organization of the 
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canids genome is the same, which indicates that the 


mitochondrial genome are evolutionary highly conserved. 


Base composition bias was determined in the red fox 
mitochondrial genome. A-T rich was found not only in 
whole genome level, but also in rRNA and tRNA genes. 
Therein, A+T frequency in tRNA genes is the highest 
one. However, the tandem repeat region determined in 
the CR is A-C rich. 

Comparison of the codon usage of red fox with 
arctic fox, gray wolf, dog and coyote suggests that they 
follow the same pattern except the unusual ATT start 
codon in the red fox. The unusual ATT start codon was 
also observed in fur seal (Lin et al, 2002) and giant 
panda (Peng, 2007). Three types of incomplete stop 
codons were determined in the red fox mitochondrial 
genome. The phenomenon of terminating with 
incomplete stop codons was also found in other 
mammals (Gissi et al, 1998; Kim et al, 1998; Ursing & 
Arnason et al, 1998; Xu et al, 1996; Peng et al, 2007; Gu 
et al, 2007). Codon bias and degeneration which occur in 
both nuclear and mitochondrial genome were found. 

The length of tRNA genes of red fox follows the 
vertebrates’ typical 59 — 75 bp. Still similar to other 
vertebrates, all tRNAs can fold into typical cloverleaf 
secondary structure except tRNASer“®”) which lacks the 
DHU stem and loop. Moreover, the anti-codons of the 22 
tRNA genes are the same as in other vertebrates. 
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